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Abstract — In  this  work,  we  analyze  the  capacity  of  airborne 
networks  where  each  node  is  equipped  with  digital  multi¬ 
beamforming  antennas.  With  this  technology,  a  node  can  form 
multiple  simultaneous  transmit  or  receive  beams  per  aperture  un¬ 
der  the  constraint  that  no  aperture  can  simultaneously  transmit 
and  receive.  Given  the  potentially  large  distances  between  air¬ 
craft,  propagation  delays  are  significant  relative  to  transmission 
times,  hence  traditional  CSMA  approaches  are  not  appropriate. 
We  first  present  our  model  of  these  multi-beamforming  capa¬ 
bilities  and  the  resulting  wireless  interference.  We  then  derive 
an  upper  bound  on  multi-access  performance  for  an  idealized 
version  of  this  physical  layer.  We  then  present  a  new  Distributed 
MAC  for  Multi-beam  Systems  (DM2S)  scheme  and  show  that 
this  random  access  protocol  achieves  the  performance  upper 
bound,  albeit  at  the  cost  of  added  delay.  We  also  consider  the 
impact  of  numerous  practical  considerations  including  stochastic 
arrivals,  latency,  and  power  constraints  on  the  performance  of 
our  random  access  MAC  protocol.  Finally,  we  present  a  system 
implementation  and  evaluation  approach  that  demonstrates  the 
feasibility  of  the  physical  layer  technology  as  well  as  the  DM2S 
scheme. 

I.  Introduction 

Today,  most  widely  deployed  airborne  communications 
waveforms  are  omni-directional  and  cannot  support  the  desired 
data  rates,  number  of  nodes,  or  the  necessary  communications 
ranges  to  support  emerging  applications  and  missions.  Further¬ 
more,  the  capacity  of  omni-directional  waveforms  is  limited 
by  the  high  amount  of  multi-user  interference.  Directional 
systems  hold  the  promise  of  higher  capacities,  increased 
scalability,  and  improved  robustness,  however,  many  direc¬ 
tional  network  implementations  fail  to  reach  the  full  potential 
afforded  by  directional  systems  due  to  their  use  of  simplistic 
networking  and  channel  sharing  schemes.  With  current  system 
designs,  the  need  for  accurate  pointing  requires  a  high  level 
of  coordination.  To  achieve  the  highest  link  data  rates,  both 
the  transmitting  node  and  receiving  node  must  simultaneously 
form  a  beam  towards  each  other.  The  difficulty  in  determining 
and  distributing  a  schedule  to  coordinate  these  transmissions 
and  receptions  in  a  mobile  ad-hoc  network  has  in  practice  led 
to  very  constrained  topologies. 

As  mentioned,  one  approach  for  system  design  is  to  de¬ 
termine  and  disseminate  a  coordination  schedule  for  node 
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transmission  and  reception  in  a  distributed  manner.  Scheduling 
in  networks  of  mobile  nodes  is  difficult  due  to  the  dynamic 
changes  in  link  states  (i.e.,  delay,  interference,  mobility),  traffic 
patterns  and  the  overhead  required  in  schedule  dissemination. 
Furthermore,  in  networks  with  large  propagation  delays,  such 
as  airborne  networks,  handshaking  based  random  access  pro¬ 
tocols  such  as  carrier  sense  multiple  access  (CSMA),  typically 
used  in  ground  networks,  are  ineffective  as  propagation  delays 
can  be  significantly  greater  than  the  packet  transmission  time. 
In  an  airborne  network,  node  separations  can  range  from  10 
-  500km,  hence  propagation  delays  can  range  from  33/zs  to  2 
ms.  We  expect  link  data  rates  of  up  to  1  Gbps  and  hence,  the 
transmission  time  of  1Kb  packet  is  on  the  order  of  l/is. 

Rather  than  carefully  coordinating  the  pointing  of  a  single 
transmit  and  receive  beam  between  pairs  of  nodes,  we  propose 
to  use  multiple  apertures  consisting  of  multi-element  digital 
antenna  arrays  and  digital  beamforming  technologies  to  form 
a  cluster  of  receive  beams  that  is  pointing  in  all  directions 
simultaneously,  thus  providing  a  continuous  360  degrees  field 
of  view  coverage  for  each  receiving  node.  This  enables  the  use 
of  low-complexity  channel  access  techniques  previously  not 
feasible  for  directional  networks,  thereby  greatly  increasing 
the  network  topology  flexibility  while  reducing  coordination 
overhead  and  hence  providing  a  practical  way  of  achieving  the 
significant  capacity  gains  possible  with  directional  networking. 
We  also  consider  employing  multiple  transmit  beams  to  pro¬ 
vide  further  gains  in  network  capacity.  These  multiple  transmit 
beams  can  be  achieved  via  multiple  array  apertures  with  one  or 
multiple  transmit  beams  per  aperture.  Multi-element  antennas 
and  digital  signal  processing  technology  have  advanced  to 
the  level  where  multiple  transmit  and  receive  beams  can  be 
formed  extremely  rapidly,  within  microseconds.  The  requisite 
neighbor  discovery  and  asynchronous  multiple  packet  acqui¬ 
sition  algorithms  that  take  advantage  of  these  capabilities 
(rapid  switching,  multiple  receive  beams,  multiple  transmit 
beams)  are  also  being  developed.  In  this  paper  we  develop 
topology  management  and  channel  access  technologies  that 
leverage  digital  arrays  and  modem  signal  processing  to  form 
a  directional  mesh  networking  system  that  provides  higher 
throughputs  and  network  capacity. 

Digital  beamforming  utilizes  an  array  of  omnidirectional 
antennas  each  equipped  with  and  Analog-to-Digital  converter 
and  with  digital  signal  processing  to  receive  and  transmit  in  an 
adaptive,  spatially  sensitive  manner.  A  digital  beam  is  formed 


Fig.  1:  Example  of  Digital  Multi-Beamforming  (DMB):  (a)  mul¬ 
tiple  simultaneous  transmit  and  (b)  multiple  simultaneous  receive 
beamforming  with  a  minimum  spatial  separation  0  —  10  degrees. 
The  center  node  cannot  form  a  transmit  beam  or  receive  beam 
toward  neighbor  A  simultaneously  with  neighbor  B  as  their  angle 
of  separation  is  within  6  degrees  of  each  other. 

on  the  receiving  end  by  applying  weights  to  the  received 
signal  vector  at  each  of  the  elements  in  the  antenna  array  to 
generate  constructive  interference  in  an  adaptive  and  spatially 
sensitive  manner  [1].  Multiple  receive  beams  in  different 
directions  can  be  formed  simultaneously.  We  assume  the 
receiver  has  the  associated  packet  acquisition  capabilities  to 
receive  and  capture  multiple  packets  from  multiple  directions 
simultaneously.  To  transmit,  a  complementary  process  is  used 
to  form  a  beam  in  the  direction  of  the  destination  [2],  [3]. 
There  are  two  mechanisms  that  enable  multiple  simultaneous 
transmit  beams.  Some  aircraft  have  multiple  apertures,  each 
with  large  multi-element  arrays.  Each  aperture  is  assumed  to 
be  able  to  independently  form  a  single  transmit  beam  or  form 
multiple  receive  beams.  Alternatively,  with  additional  signal 
processing,  a  single  aperture  can  be  used  to  simultaneously 
form  beams  to  multiple  receive  nodes.  Each  beam  potentially 
carries  independent  information.  The  total  transmit  power 
is  fixed,  hence  forming  multiple  transmit  beams  results  in 
splitting  the  power  between  the  multiple  beams.  Given  these 
assumptions  for  the  physical  layer  beamforming  capabilities, 
we  develop  and  assess  technologies  for  sharing  the  wireless 
channel.  We  begin  by  calculating  and  upper  bound  on  the 
capacity  of  a  network  with  these  multiple  beam  capabilities. 
We  then  design  a  random  access  protocol  and  show  that  this 
protocol  achieves  network  capacity. 

II.  Related  Work 

There  have  been  numerous  studies  on  the  throughput  of 
single  beam  directional  networks  [1],  [4]— [7].  However,  many 
of  these  studies  are  protocol-based  and  evaluated  via  discrete 
event  simulators;  little  analytical  work  has  been  done  for  single 
beam  directional  networks  [8]— [10] . 

There  are  also  works  which  analyze  the  throughput  of 
networks  with  multi-packet  reception  (MPR)  [11],  [12]  using 
802.11  WiFi  [13],  Code  Division  Multiple  Access  (CDMA) 
[14],  Multiple-Input  Multiple-Output  (MIMO)1  [16].  However, 
they  do  not  consider  multi-packet  transmission  as  well. 

Finally,  existing  studies  which  consider  multi-transmit  and 
multi-receive  beamforming  capabilities  are  access-point  based 

Unlike  beamforming,  MIMO  utilizes  multiple  antennas  to  form  numerous 
low  data  rate  signals.  Additionally,  MIMO  works  focus  on  interference 
cancellation  (e.g.,  [15]),  and  are  out  of  the  scope  of  this  work. 


and  utilize  carrier-sense  multiple  access  (CSMA)  or  RTS/CTS- 
based  schemes  [3],  [17]— [19]  or  require  coordination  (i.e.,  two- 
hop  neighbors)  [20]-[22]. 

A  review  of  adaptive  antenna  array  techniques  which  can 
perform  multi-packet  transmission  or  reception  can  be  found 
in  [2].  However,  following  a  common  assumption,  the  beam¬ 
forming  techniques  are  assumed  to  be  too  complex  for  mobile 
nodes  and  are  reserved  to  the  basestation.  Unlike  previous 
work,  in  this  work,  we  consider  networks  with  (1)  large 
propagation  delays  between  nodes  making  RTS/CTS  schemes 
ineffective  and  (2)  homogenous  devices  with  multi-beam  ca¬ 
pabilities  as  outlined  in  Section  III-A. 

III.  Model  and  Problem  Definition 

In  this  section  we  outline  the  network  model,  assumptions, 
and  performance  metrics  as  well  as  the  throughput  problem 
definition. 

A.  Beamforming  Capabilities 

In  this  work,  each  node  is  equipped  with  a  digital  beam¬ 
forming  antenna.  There  are  numerous  benefits  to  using  beam¬ 
forming  antennas  (e.g.,  improved  transmission  reliability  and 
transmission  range  extension  [1]).  However,  in  this  work,  we 
focus  on  analyzing  the  capacity  improvements  stemming  from 
the  increased  spatial  reuse  of  forming  multiple  beams.  As  such 
we  will  leverage  the  following  capabilities: 

Observation  1:  [Multi-beam  transmit,  Fig.  1(a)]  A  node  can 
transmit  up  to  F?max  simultaneous  directional  beams  of  width 
6  degrees  to  their  neighbors. 

As  described  in  the  section  above,  a  node  can  radiate  power 
in  multiple  directions,  forming  distinct  simulatenous  transmis¬ 
sions.  However,  this  process  is  limited  by  the  angle  of  the 
radiated  power,  termed  the  beamwidth ,  and  denoted  by  6.  As 
such,  a  node  cannot  form  two  transmissions  within  0  degrees 
of  one  another.  Additionally,  this  process  is  limited  by  the 
processing  and  power  constraint  of  the  node.  Hence,  each  node 
is  limited  to  forming  Bmax  simultaneous  transmission  beams. 

Observation  2:  [Multi-beam  receive,  Fig.  1(b)]  A  node  can 
receive  simultaneous  directional  beams,  provided  they  are 
separated  by  at  least  0  degrees. 

On  the  receiving  node,  a  transmission  can  only  be  discerned 
from  that  of  a  neighbor  if  they  are  separated  by  6  degrees. 
In  the  case  of  two  or  more  beams  arriving  within  0  degrees, 
we  say  that  none  of  the  beams  can  be  correctly  decoded. 
It  is  important  to  note  that,  using  digital  signal  processing 
techniques  to  constantly  sample  and  decode  the  antenna  array, 
this  allows  digital  beams  to  be  formed  on  the  receiving  end 
of  the  antenna  after  a  packet  has  already  arrived  at  a  node. 
This  effectively  enables  a  digital  beamforming  antenna  to  act 
in  an  omnidirectional  mode,  but  have  antenna  gains  equivalent 
to  directional  mode. 

Observation  3:  [No  transmit  while  receive]  A  node  cannot 
simultaneously  form  a  transmit  beam  while  forming  a  receive 
beam. 

This  third  and  final  observation  stems  from  the  fact  that  a 
transmitting  beam  will  interfere  with  any  receiving  beam  in  the 
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antenna  antenna  array.  We  will  assume  that  anytime  a  message 
arrives  while  a  node  is  transmitting  a  beam,  that  message  is 
lost.  However,  it  should  be  noted  that  the  message  which  is 
transmitted  can  still  be  correctly  decoded  at  it’s  destination. 

B.  Network  Model 

There  are  n  mobile  airborne  nodes  (assumed  to  be  on  a 
2  dimensional  plane)  and  all  nodes  are  within  range  of  each 
other.  However,  as  beamforming  nodes  can  be  quite  far  from 
one  another  and  still  be  within  communication  range2,  there 
can  be  significant  propagation  delays. 

We  assume  that  all  nodes  have  pointing,  acquisition,  and 
tracking  knowledge  of  each  others  whereabouts  such  that  they 
can  determine  the  angle  of  a  beam  to  a  neighbor.  Due  to  the 
distance  between  the  nodes  and  a  small  transmit  angle  (e.g., 
10°),  even  relative  information  of  a  neighbors  location  is  easy 
to  point  and  track. 

Each  node  is  equipped  with  a  DMB  antenna  (see  Sec¬ 
tion  III- A).  We  will  assume  a  packet  network,  with  fixed 
unit  length  transmission  times.  Nodes  can  have  two  types  of 
traffic:  backlogged  or  stochastic  arrivals.  Initially,  to  analyze 
the  capacity  of  the  network,  we  will  first  consider  the  case 
where  all  nodes  have  an  infinite  buffer  of  packets  to  send  to 
each  of  their  neighbors.  In  Section  V-A,  we  will  relax  this 
assumption  and  consider  nodes  which  have  packet  arrivals 
according  to  a  poisson  process  with  arrival  rate  A. 

C.  Problem  Definition 

A  node  wishes  to  maximize  the  traffic  being  sent  to  each 
of  their  n  —  1  neighbors  (all-to-all  traffic).  As  a  performance 
metric,  we  will  consider  the  per-link  throughput. 

Definition  1:  For  a  directed  link  between  a  source  and  a 
destination,  the  per-link  throughput  is  the  percentage  of  time 
that  the  source  is  transmitting  and  the  packet  is  successfully 
received  at  the  destination. 

Note  that  as  this  is  a  fully  connected  graph,  there  are  n2  —  n 
directed  links  in  the  network. 

Fundamentally,  this  problem  is  made  difficult  by  the  lack 
of  coordination  causing  each  node  to  act  independently.  The 
throughput  is  composed  of  the  rate  that  a  node  transmits  (pack¬ 
ets/second)  and  multiplied  by  the  likelihood  that  the  packet 
is  successfully  received  at  the  destination.  To  be  successfully 
received,  the  destination  node  must  not  be  transmitting  any 
packets  for  the  duration  of  the  packet  arrival.  As  the  packet 
will  arrive  at  the  destination  after  some  unknown  (and  po¬ 
tentially  long)  propagation  delay,  it  is  very  difficult  to  have 
the  destination  node  anticipate  a  packet  arrival.  Therefore,  by 
having  a  node  transmit  often,  it  will  result  in  lower  success 
rates.  Hence,  there  is  a  tradeoff  for  attempting  to  transmit 
often  with  numerous  collisions  or  transmitting  infrequently 
with  fewer  collisions.  In  the  next  section,  we  will  demonstrate 
an  uncoordinated  stochastic  protocol  which  can  achieve  the 
optimal  network  throughput,  even  if  it  were  coordinated. 

2The  range  of  airborne  beamforming  nodes  can  be  up  to  500km,  as  all 
nodes  have  Line-of-sight  (LOS)  to  each  other. 


n 

Number  of  nodes  within  communication  range 
of  one  another 

0 

Beam  width  size  (degrees) 

Bm&x 

Maximum  number  of  transmit  beams 

T 

The  per-link  throughput 

IV.  Throughput  Analysis 

We  now  begin  the  throughput  capacity  analysis  for  a  net¬ 
work  of  n  nodes  with  DMB  antennas  whose  capabilities  were 
outlined  in  Section  III-A.  To  understand  the  upperbounds  on 
throughput  capacity,  in  this  section  we  consider  an  idealized 
physical  layer  with  no  beam  width  or  power  constraints. 

In  the  analysis  below,  we  first  derive  an  upperbound  on  the 
capacity  of  such  a  system.  Following  this,  we  show  that  a 
random  access  MAC  can  also  achieve  the  same  capacity  as 
the  scheduled  MAC  in  the  asymptotic  limit  of  delay. 

A.  Throughput  Upperbound 

As  an  upperbound,  we  will  consider  the  throughput  of 
the  muti-beam  sectorized  directional  system  under  relaxed 
versions  of  Observations  (1)  and  (2)  that  nodes  can  form 
any  number  of  beams  (F>max  =  n  —  1)  of  arbitrarily  small 
angle  ( 0  =  0).  Therefore,  the  main  constraint  stems  from 
Observation  (3)  that  a  node  cannot  simultaneously  transmit 
and  receive  a  packet.  To  derive  the  upperbound,  we  will 
assume  that  the  nodes  are  perfectly  synchronized  and  have 
aligned  slots. 

Intuitively,  when  considering  scheduling  in  such  a  system, 
the  problem  reduces  to  a  graph  vertex  coloring.  Put  simply, 
if  there  exists  a  vertex  coloring  with  only  2  colors  (i.e.,  the 
graph  is  bi-partite),  then  the  optimal  scheduling  is  to  schedule 
all  of  the  black  nodes  in  a  single  slot  to  transmit  to  all  of  their 
neighbors,  each  of  which  are  red  nodes.  In  the  subsequent 
slot,  schedule  all  of  the  red  nodes.  In  this  scenario,  clearly 
each  link  will  be  activated  every  other  slot,  and  all  packets 
will  be  successfully  received.  Hence,  the  throughput  per  link 
is  50%  and  this  is  an  upperbound.  However,  it  is  clear  that 
such  an  upperbound  will  not  be  able  to  be  achieved  in  practical 
scenarios  where  network  topology  is  not  a  bi-partite  graph. 

For  a  complete  graph  (as  described  in  Section  III),  we  can 
also  determine  an  upperbound  on  throughput. 

Theorem  1:  The  upperbound  on  per-link  throughput  is, 

T*  =  25%.  (1) 

Proof:  With  nodes  synchronized,  slotted,  and  no  propa¬ 
gation  delay,  the  throughput  in  a  slot  is  completely  determined 
by  the  set  of  nodes  which  are  transmitting  vs.  the  set  of 
nodes  which  are  receiving.  As  a  simple  example,  if  one  node 
transmits  in  a  slot  to  the  remaining  n  —  1  nodes,  then  there 
are  n  —  1  successfully  activated  links  in  that  slot.  Following 
this,  if  j  nodes  transmit  in  a  slot  to  n  —  j  nodes  which  receive, 
the  number  of  successfully  activated  links  will  be  j  *  (n  —  j ) . 
It  is  easy  to  see  that  an  upperbound  on  throughput  will  then 
be  when  n/2  nodes  transmit  to  n/2  nodes  for  an  aggregate  of 


Distributed  MAC  for  Multibeam  Systems  (DM2S) 

Input:  Beamwidth  0 ,  Power  constraint  Bmax,  Node  i ,  Angles  ay  to 
all  neighbors  j 

1:  while  Queue  of  waiting  packets  is  not  empty  do 
2:  5  =  0 

3:  for  each  neighbor  j  do 

4:  if  | Oii  —  a.j |  >  G  5  and  |5|  <  Bmax  then  S  =  S  U  i 

5:  Form  a  transmit  beam  to  all  neighbors  in  S. 

6:  Enter  into  receive  mode  for  a  random  period  of  time  R 
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Fig.  2:  Random  Access  MAC  Overview:  After  waiting  a  random 
period  of  time,  R ,  a  node  transmits  a  unit-length  packet  on  it’s  link. 
If  the  entire  packet  arrives  within  the  stochastic  interarrival  time  R , 
it  is  successfully  received. 


n2  / 4  activated  links.  As  there  are  «  n2  links,  this  results  in 
a  per  link  throughput  upperbound  of  2 5%. 3  ■ 

Again,  however,  this  scheme  requires  that  nodes  are  per¬ 
fectly  synchronized  (i.e.,  aligned  slots)  and  coordinated  which 
is  difficult,  if  not  impossible,  to  achieve  in  practice.  Nonethe¬ 
less,  in  the  following  section,  we  will  propose  and  analyze  a 
number  of  schemes  which  will  approach  this  upperbound  of 
25%  throughput  per  link. 

In  the  section  above,  we  demonstrated  that  25%  activa¬ 
tion  per  link  is  an  upperbound  on  the  throughput.  In  the 
networks  of  interest,  a  scheduled  system  does  not  work  as 
coordination  is  too  difficult  to  achieve  in  an  ad  hoc  network. 
This  is  amplified  in  airborne  networks  where  propagation 
delays  are  non-negligible.  Furthermore,  carrier  sense  multiple 
access  (CSMA)  is  not  practical  as  propagation  delays  cause 
inaccurate  channel  assessment  and  the  RTS/CTS  message 
exchange  becomes  prohibitive.  Therefore,  in  this  subsection, 
we  will  demonstrate  a  class  of  protocols  which  can  achieve 
this  without  any  coordination  between  the  nodes. 

B.  Distributed  MAC  for  Multibeam  Systems  (DM2S) 

In  this  section,  we  now  describe  an  Aloha-esque  MAC  layer 
protocol  which  doesn’t  require  any  coordination  between  the 
nodes. 

To  send  a  packet,  a  node  will  wait  a  random  amount  of  time 
R ,  then  enter  transmitting  mode.  While  transmitting,  the  node 

3  A  similar  problem  was  considered  [22].  However,  they  use  a  tree  structure 
to  select  the  set  of  nodes  to  transmit  and  receieve  in  each  slot,  providing 
guarantees  of  link  activations  in  a  superframe.  For  completeness,  we  have 
recreated  a  similar  argument. 


will  send  a  packet  to  all  of  it’s  neighbors  that  it  can  (depending 
on  n  and  T?max).4  During  the  waiting  time,  a  node  is  in  receive 
mode  and  can  successfully  receive  a  packet  from  any  of  it’s 
neighbors  if  it  arrives  entirely  within  the  receive  window  (i.e., 
it  does  not  overlap  with  the  start  or  end  of  the  interval).  After 
transmitting  the  packet,  the  node  again  waits  a  random  amount 
of  time  before  sending  the  next  packet.  The  probability  density 
function  and  cumulative  distribution  function  are  represented 
as  /i?(r)  and  FR(r),  respectively. 

This  abstraction  of  a  MAC  protocol  contains  many  vari¬ 
ations  of  existing  MAC  protocols.  For  example,  the  well 
known  slotted  Aloha  can  be  thought  of  as  having  R  distributed 
according  to  a  geometric  distribution  with  discrete  unit-length 
intervals  and  probability  p  transmitting  in  a  given  slot. 

For  the  scheme  where  the  distribution  of  the  receive  window 
time  is  arbitrary.  We  now  derive  the  throughput  per  link  in  this 
scenario: 


^dm2s  — 


E[R] 


1+E[i?]  1  +  E  [R] 


P (Y(t)  >  1). 


(2) 


The  first  term  represents  the  average  time  that  a  link  is 
transmitting.  The  2nd  term  represents  the  average  time  the 
receiver  is  in  receive  mode.  In  the  3rd  term,  Y  (t)  represents 
the  forward  recurrence  time  of  the  receive  slot  given  an  arrival 
at  time  t.  The  3rd  term  represents  the  probability  that,  given 
a  packet  arrives  while  in  a  receive  slot,  the  receive  slot  has  a 
remaining  length  greater  than  the  length  of  the  packet  being 
transmit. 

Using  well  known  results  from  probability  theory  [23],  we 
have 


E[F(t)]  = 


E[i?2] 

2E[i?] 


E[R2} 
2  L 


(3) 


l 

Wt)  >  i)  =  1  -  j 1  -  FR(r)dr  (4) 
0 

Equation  (3)  shows  that  the  expected  length  of  the  forward 
recurrance  time  grows  with  the  second  moment  of  the  receive 
window  size.  Equation  (4),  which  is  valid  only  for  non-lattice 
distributions  suggests  that  for  F(Y(t)  >  L)  -A  1,  FR{r)  «  1 
for  some  small  value  of  r. 

Using  the  insights  above,  in  the  following  sections,  we  con¬ 
sider  a  couple  of  variations  on  the  distribution  of  the  receive  in¬ 
terval  duration  and  quantify  the  corresponding  throughput.  Our 
goal  is  to  find  distributions  which  maximize  the  throughput. 
Note  that,  by  ignoring  the  3rd  term  in  (2)  and  attempting  to 
maximize  the  first  two  terms,  we  set  E  [R\  =  1.  This  means  that 
we  would  be  transmitting  for  50%  of  the  time  and  receiving 
for  50%  of  the  time  for  a  throughput  of  25%,  which  is  the 
upperbound.  However,  to  maximize  the  3rd  term  such  that 
P( Y(t )  >  1),  using  (4),  we  must  show  that, 


l 


/ 


FR(r)dr  ->■  1. 


transmissions  cannot  cause  collisions  with  other  links.  Hence,  if  a  node 
is  transmitting,  it  will  transmit  to  as  many  neighbors  as  possible. 


Intuitively,  this  means  that  despite  having  an  expected  value 
at  1,  nearly  all  of  the  weight  of  the  distribution  of  R  must  fall 
before  1.  This  observation  will  be  leveraged  later  on. 

Before  showing  a  formal  proof  characterizing  the  class  of 
capacity -achieving  algorithms,  we  will  compute  the  through¬ 
put  for  commonly  used  distributions. 

Exponential  Distributed:  In  this  case,  let  the  receive  window 
length  R  ~  exponential  with  mean  E[R].  Then  the  throughput 
can  be  written  as 


T  2  _  L  E[R]  —L/E[R] 

Jdm2s-exp  -  E[R]+LE[R]+Le 

where  E  [R\  is  the  expected  duration  of  a  receive  window.  The 
first  term  represents  the  average  time  that  a  link  is  transmitting. 
The  2nd  term  represents  the  average  time  the  receiver  is  in 
receive  mode.  The  3rd  term  represents  the  probability  that, 
given  a  packet  arrives  while  in  a  receive  slot,  the  receive  slot 
has  a  length  greater  remaining  than  the  length  of  the  transmit 
slot.  Note  that,  for  the  exponential  distribution,  the  memoryless 
property  ensures  that  the  distribution  of  the  remaining  time  in 
an  interval  is  equivalent  to  the  distribution  of  the  size  of  the 
interval. 

This  value  is  maximized  at  E\R]  =a? 1  +  y/2  for  a  value  of 
throughput, 

Ce_1/c 

^dm2s-exp*  =  ^  (jy  =  13.68% 
where  C  —  1  +  \[2. 

Slotted  and  Unsynchronized  In  this  scheme,  the  interarrival 
time  of  packets  is  selected  according  to  gemoetric  distribution 
with  discrete  intervals  of  length  1.  That  is, 


p(l-p)r  :Vr  =  0,1, 2, 3... 
0  :  else 


In  this  scheme,  all  nodes  divide  time  into  slots  where  the 
slot  length  is  the  duration  of  time  required  to  send  a  packet. 
In  each  slot,  a  node  will  transmit  with  probability  p.  Hence, 
slotted  aloha-like  analysis  can  be  utilized  to  determine  the 
throughput,  T. 


Tdm2s-slot  =  ^(Src  Tx)P(Dest  Rx  for  2  Slots)  (5) 
=  P{  1  -  P )2  (6) 

Clearly,  the  throughput  is  maximized  at  p  =  1  /3  for  a  value 
of 

^dm2s-slot*  —  ^  —  14.81%  (7) 

C.  Capacity  Achieving  Distribution 

In  the  above  scenarios,  we  demonstrated  that  the  expo¬ 
nential  and  geometric  distributions,  when  optimized,  result  in 
a  per  link  throughput  of  13.68%  and  14.81%,  respectively. 
We  hypothesize  that  the  added  performance  for  the  slotted 
distribution  over  the  exponential  distribution  stems  from  the 
added  variance  of  the  length  of  the  receive  window  in  the 


latter  case.  As  alluded  to  in  Section  IV-B,  the  2nd  moment 
of  the  receive  window  is  a  main  factor  in  the  likelihood  that 
a  node  will  interrupt  an  incoming  transmission  by  beginning 
the  transmit  mode.  Hence,  in  this  section  we  consider  a 
MAC  scheme  which  leverages  large  variances  in  the  receive 
window  size.  A  similar  idea  was  initially  proposed  in  [24], 
albeit  for  omnidirectional  networks  and  no  analysis  of  capacity 
achieving  distributions  is  given. 

Theorem  2:  There  exists  a  distribution  for  the  receive  win¬ 
dow  size  which  results  in  the  random  access  MAC  scheme 
asymptotically  achieving  the  upperbound  on  per-link  through¬ 
put. 

Proof:  We  will  prove  this  theorem  through  construction 
of  the  Capacity  Achieving  Distribution  (CAD)  which  has  a 
receive  window  size  distribution  as  follows: 


:  r  =  0 
:  r  =  x  ’ 


where  x  is  referred  to  as  the  backoff  time  parameter  of  the 
distribution  with  the  constraint  that  x  >  1.  Recall  that  1  is 
the  duration  of  a  packet  transmission  time.  In  this  distribution 
there  is  a  (relatively)  large  probability  that  the  receive  window 
size  is  of  length  0,  imply  back-to-back  packet  transmissions. 
On  the  other  hand,  there  is  a  relatively  small  probability  that 
the  receive  window  is  of  length  x ,  where  x  is  arbitrary.  Hence, 
we  will  refer  to  x  as  the  delay  parameter  as  it  will  represent 
the  amount  of  time  a  node  must  wait  for  a  receive  interval  to 
end  to  send  a  packet. 

Therefore,  in  this  distribution,  the  mean  receive  window 
length  is  E  [R\  =  1  and  the  average  time  receiving  is  equal 
to  the  average  time  transmitting.  The  second  moment  of  the 
receive  window  is  E  [R2]  =  x.  Therefore,  in  this  distribution, 
the  mean  independent  of  x,  but  the  variance  clearly  can  grow 
to  infinity. 

For  a  given  value  of  x,  the  probability  that  the  forward 
recurrence  time  is  greater  than  the  length  of  the  packet  at  an 
arbitrary  time  t  is 

p(y(t)  >  i)  =  f 

X 

Therefore, 

1  i  x  —  1 

lim  ^~dm2s-cad  =  lim  —  —  =  25%, 

x^-oo  UM  5  x — >-00  2  2  X 

and  the  overall  throughput  of  the  system  approaches  25%.  The 
relationship  between  the  throughput  and  the  delay  parameter, 
x ,  is  shown  in  Fig.  3(a).  ■ 

Following  this  theorem,  we  now  show  that  infinite  variance 
is  a  requirement  for  any  capacity-achieving  random  access 
MAC  scheme. 

Lemma  1:  For  any  receive  window  size  distribution  with 
finite  variance  and  unit  mean,  the  DM2S  scheme  does  not 
achieve  the  capacity  of  25%. 

Proof:  Clearly,  from  (2)  the  throughput  can  only  be  25% 
when  P (Y(t)  >  1)  =  1.  For  any  finite  variance  distribution, 
from  (4),  it  is  clear  that  P (Y(t)  >  1)  <  1  and  therefore  the 
25%  capacity  cannot  be  achieved.  ■ 
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Gamma  Shape  Parameter 
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Pareto  Shape  Parameter 

(b) 


Fig.  3:  a)  The  throughput  for  the  capacity  achieving  distribution 
approaches  the  capacity  of  25%  as  the  delay  parameter  x  — ►  oo. 
b)  Random  access  MAC  per-link  throughput  with  R  distributed 
according  to  the  Pareto  and  Gamma  distributions. 


D.  Infinite  Variance  Distributions 

As  demonstrated  in  the  above  section,  a  receive  window 
distribution  with  finite  mean  and  infinite  variance  will  allow 
the  MAC  scheme  to  approach  the  upperbound  capacity  of 
25%  per-link  throughput.  In  this  subsection,  we  show  exper¬ 
imentally  that  this  holds  for  two  well  known  distributions: 
the  Pareto  and  the  gamma  distributions.  In  Fig.  3(b),  we 
select  the  scale  parameters  for  both  distributions  such  that  the 
mean  receive  window  length  is  equal  to  1.  Then,  we  vary  the 
shape  paraemters.  That  is,  for  the  pareto  distribution  the  scale 
parameter  =  (a  —  1  )/a  where  a  is  the  shape  parameter. 
For  the  gamma  distribution,  the  scale  parameter  is  1  /k  where 
k  is  the  shape  parameter.  As  shown  in  the  figure,  as  we  vary 
the  shape  parameters  towards  0,  the  corresponding  throughput 
approaches  the  capacity  of  25%. 

For  the  gamma  distribution,  it  is  expected  that  as  the 
shape  parameter  approaches  0,  the  variance  of  the  distribution 


approaches  infinity  (importantly,  with  finite  mean).  This  is 
consistent  with  the  conjecture  that  it  is  indeed  the  variance 
which  is  key  towards  approaching  the  capacity.  However,  in 
the  pareto  distribution  with  the  shape  parameter  1  <  a  <  2, 
the  receive  window  size  has  infinite  variance,  independent  of 
the  shape  parameter.  Hence,  we  find  that  despite  having  a 
distribution  with  infinite  variance,  we  do  not  necessarily  have 
a  capacity  achieving  distribution.  Instead,  the  throughput  of 
the  random  access  scheme  utilizing  the  pareto  distribution  ap¬ 
proaches  the  25%  capacity  as  the  shape  parameter  approaches 
0.  We  hypothesize  that  is  not  the  variance,  but  the  2nd  moment 
growing  to  infinity  which  is  sufficient  condition  for  a  capacity 
acheiving  MAC  scheme.  Unfortunately,  a  proof  of  such  a 
conjecture  has  thusfar  been  out  of  reach  and  is  a  subject  of 
future  work. 

V.  Practical  Considerations 

In  the  prior  Section,  we  analyzed  and  demonstrated  that 
the  upperbound  per-link  throughput  could  be  achieved  via  a 
random  access  MAC  protocol  under  idealized  physical  layer 
assumptions  where  there  was  an  infinite  queue  of  packets,  no 
beamwidth  ( 0  =  0),  and  no  power  constraints  (Bmax  =  n  —  1). 
In  this  section,  we  relax  those  assumptions  in  turn. 

First,  we  consider  the  scenario  where  packets  arrive  accord¬ 
ing  to  a  poisson  arrival  process.  Then,  we  consider  the  case 
when  the  the  spatial  separation  between  nodes  must  be  greater 
than  0  to  transmit  or  to  receive  from  a  neighbor.  Last,  we 
consider  the  per-link  throughput  in  the  presence  of  a  power 
constraint. 

A.  Stochastic  Arrivals 

In  this  subsection,  we  relax  the  constraint  of  all  nodes 
having  an  infinite  buffer  of  packets  to  send.  Instead,  packets 
arrive  according  to  a  stochastic  process  with  arrival  rate  A 
to  each  node.  When  a  packet  arrives,  a  distinct  packet  is 
generated  for  each  of  the  n  —  1  neighbors  of  a  node. 

This  system  can  then  be  analyzed  as  an  M/G/l  queue  [25] 
which  has  a  steady-state  probability  of  having  a  non-empty 
queue  of  ttq  =  1  —  A / fi  where  A//x  is  the  utilization  of  the 
queue  based  on  an  arrival  rate  A  and  average  service  rate  p. 

To  determine  the  throughput,  we  will  utilize  a  similar 
approach  as  Section  IV.  First,  we  find  the  percentage  of  time 
that  a  given  link  is  activated  by  the  source.  Then,  we  find  the 
probability  that  a  given  packet  is  successfully  received.  The 
MAC  scheme  is  as  per  the  above  section. 

A  node  will  transmit  if  it  has  at  least  1  packet  in  the  queue 
and  it  is  in  a  transmission  window.  Hence,  it  transmits  with 
rate  (1  —  tto)*1/(1A-E[R])  =  p*  1/2.  A  packet  is  successfully 
received  if  the  receiving  node  a)  has  an  empty  queue  with  no 
new  arrivals  for  a  packet  duration,  b)  has  an  empty  queue 
with  an  arrival  sometime  in  the  next  unit  time  but  decides  to 
wait  to  transmit,  or  c)  has  a  non-empty  queue  and  is  in  the 
receive  mode  for  length  at  least  unit  time.  These  3  events  are 
independent  and  the  probability  of  success  can  be  written  as: 

^ r0e_A  +  7r0(l  -  e-A)  *  l/x  +  (1  -  7r0)  *  - - . 

2  x 


Then  the  overall  throughput  is 
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^dm2s-cad 


p2  x  —  1  p(  1  —  p) ,  X 
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1  —  e  A 
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(8) 


It  is  easy  to  see  that  as  p  — >>  l,x  — )>  oo,  the  queues 
become  backlogged  and  the  throughput  approaches  the  same 
25%  upperbound  throughput  from  Section  IV- A.  However,  this 
model  allows  us  to  observe  the  tradeoff  for  a  stochastic  arrival 
process  between  throughput  and  delay. 

To  compute  the  delay,  we  will  utilize  the  well-known 
Pollaczek-Khinchin  formula  [25]  which  states  that  the  total 
delay  in  system  is 


Delay  =  Y  + 


AY2 


2(1  -PV 
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where  Y  and  Y2  is  the  average  service  time  and  second 
moment  of  service  time,  respectively.  The  distribution  of 
service  time  of  a  packet  is: 


fy(y) 


■  y  =  i 

:  y  =  x  +  1 


Recall  that  x  is  the  parameter  of  the  capacity  achieving 
distribution.  The  service  time  distribution  is  simply  the  delay 
from  the  MAC  protocol  with  the  added  unit-time  required 
to  transmit  the  packet.  Then,  E [Y]  =  =  2  and 

E[Y2}  =  =  3  +  x.  The  delay  can  then  be  written 

as: 


Delay  =  2  + 


A(x  H-  3) 
2(1  -  A/2)  * 


(9) 


Fig.  4  shows  the  relationship  between  throughput  (8)  and  delay 

(9). 

Fig.  4(a)  shows  the  throughput  for  various  utilization  rates 
p.  For  larger  values  of  p ,  the  throughput  improves  with  larger 
backoff  times  as  the  system  behaves  similar  to  a  system  with 
a  queue  that  never  empties.  Corresponding  to  a  system  with 
an  infinte  queue  of  packets,  as  p  1  and  x  — >■  oo,  the  system 
approaches  the  upperbound  throughput  of  25%.  However,  for 
smaller  values  of  p ,  the  best  scheme  is  to  simply  have  no 
backoff  time  and  to  immediately  transmit  a  packet  when  it 
arrives. 

In  Fig.  4(b),  p  is  varied  for  various  values  of  x  and  the 
corresponding  delay  D  from  (9)  and  throughput  from  (8)  are 
shown.  Again,  for  smaller  values  delay,  we  see  that  larger 
values  of  backoff  time  result  in  larger  throughput.  However, 
at  some  threshold  it  becomes  beneficial  to  increase  the  backoff 
time  parameter.  Also,  note  that  the  delay  increases  linearly  in 
a  logarithmic  scale  with  the  throughput  for  small  values  of 
delay.  There  is  an  elbow  in  the  curve  around  103  slots  at 
which  point  the  growth  to  25%  throughput  is  relatively  slow. 
This  implies  that  a  practical  scheme  may  settle  for  slightly 
less  than  optimal  throughput  for  better  delay  performance 


Fig.  4:  Relationship  between  per-link  throughput  (TDM2S_CAD)  and 
delay  ( D ):  a)  Link  throughput  as  a  function  of  the  backoff  time 
parameter  of  the  CAD  for  the  DM2S  scheme  and  each  node  has  a 
utilization  of  p,  b)  Link  throughput  as  a  function  of  the  delay  (in 
slots)  for  various  values  of  the  backoff  time  parameter  x. 


B.  Beamwidth  (0) 

From  the  perspective  a  receiving  node,  neighbors  must  be 
spatially  distributed  by  at  least  0  degrees  in  order  to  decode 
the  separate  streams.  Thereby,  we  now  condition  on  the  case 
where  a  destination  node  is  receiving  from  some  source  node 
according  to  the  capacity  achieving  MAC  protocol  described 
above.  To  facilitate  this  analysis,  we  will  now  assume  that 
nodes  are  uniformly  distributed  according  to  a  two  dimensional 
spatial  poisson  process  with  mean  n  neighbors  per  square 
unit  area.  Thus  far,  we  have  not  been  able  to  quantify  what 
the  upperbound  on  per-link  throughput  would  be  given  the 
beamdwidth  constraint.  However,  we  will  now  derive  the 
expected  per  link  throughput  for  the  DM2S  protocol. 

In  this  case,  we  focus  our  analysis  on  a  given  link  between 
to  a  neighbor.  The  stream  can  correctly  be  decoded  if  no  other 
neighbors  within  0  degrees  is  transmitting  for  duration  of  a 


Fig.  5:  Relationship  between  per- link  throughput  ( T9 )  and  node 
density  n  when  nodes  are  distributed  according  to  a  spatial  poisson 
process  and  have  a  beamdwidth  constraint  0. 


packet.  To  determine  the  expected  drop  in  throughput,  we 
define  the  set  of  number  of  neighbors  within  6  degrees  as 
a  random  variable  we  will  denote  as  Z.  For  this  analysis  we 
will  assume  that  the  DM2S-CAD  protocol  is  being  used  with 
the  backoff  time  parameter  x  =  oo. 

-^dm2s-cad  =  Tdm2s-cadp(^  rx’ing  f°r  1  slot) 

1  OO 

=  -  rx’ing  | Z  =  z)¥(Z  =  z) 

^  z= 0 

_  1  y*,l)x(n0/36O)z  0/36O 

zl 

z=0 

_  \  —  n0/360 
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Fig.  6:  Relationship  between  per-link  throughput  and  number  of 
transmit  beams  Bmax\  a)  upperbound  per-link  throughput  (T*,Smax) 
compared  to  DM2S  per-link  throughput  (T_^T  ),  b)  approxima- 

DM  o-LAU 

tion  ratio  between  T*,Smax  and  TSrT 

dm2s-cad 


This  result  is  quantified  in  Fig.  5. 

C.  Power  (Bmax) 

In  this  subsection,  we  assume  that  there  is  a  power  constraint 
limiting  the  number  of  transmit  beams  that  a  node  can  form 
to  T>max.  Using  a  similar  argument  as  Section  IV- A,  we  derive 
the  upperbound  of  the  per  link  throughput  as 


x/n  and  the  throughput  for  the  DM2S  protocol  can  then 
be  derived  as 

TDm“s-CAD  =  TDM2S-CADBma  x/n. 

The  resulting  throughput  and  approximation  ratios  can  be 
found  in  Fig.  6.  The  approximation  ratio  varies  between  4 
(at  T>max  =  1)  to  2  (at  Bmax  =  n/2)  to  1  (at  Bmax  =  n). 


-^max(^  -^max)/^  •  -^max  ^  Tl/2 

25%  :  5max  >  n/2  ' 

As  can  be  seen,  the  upperbound  throughput  with  the  power 
constraint  is  identical  to  the  upperbound  without  a  power 
constraint  if  T>max  >  n/2.  This  is  because  in  the  upperbound 
solution,  each  transmitting  node  sends  to  n/2  of  their  neigh¬ 
bors,  hence  they  only  need  form  n/2  transmit  beams. 

For  the  DM2Sprotocol,  we  then  adapt  it  such  that  when 
transmitting  to  neighbors,  choose  the  set  of  T>max  neighbors 
to  transmit  to  randomly.  Therefore,  when  transmitting,  a 
transmission  beam  will  be  sent  to  a  neighbor  with  probability 


VI.  Conclusions 

We  consider  a  communications  network  of  airborne  plat¬ 
forms  where  each  node  is  equipped  with  digital  multi¬ 
beamforming  antennas.  With  this  technology,  a  node  can  form 
multiple  simultaneous  transmit  or  receive  beams  under  the 
constraint  of  no-transmit-while-receive.  We  first  present  an 
idealized  model  of  these  physical  layer  capabilities.  Using  this 
model,  we  derive  an  upper  bound  on  the  network  capacity 
assuming  an  optimal  MAC  protocol.  We  then  present  the 
Distributed  MAC  for  Multi-beam  Systems  (DM2S)  scheme 
and  show  that  this  random  access  protocol  achieves  the  perfor¬ 
mance  bound.  We  assess  the  resulting  delay  implications.  We 


then  relax  our  initial  idealistic  assumptions  and  consider  the 
impact  of  numerous  practical  considerations  including  power 
constraints,  beam  widths,  stochastic  arrivals,  and  latency  re¬ 
quirements  on  the  performance  of  the  new  random  access 
MAC  protocol. 
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